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Abstract. We investigated the optical properties of 7 clusters of galaxies observed in three colors over the range of absolute 
magnitudes -24 < M < -12. Our aim is to estimate the Luminosity Function and the total cluster luminosity of our sample 
in order to have information about the formation and the evolution of galaxies in clusters. In this paper, we present the main 
points of our analysis and give the formal parameters obtained by fitting the data using the maximum likelihood algorithm. We 
find consistency between our results and other works in literature confirming the bimodal nature of the luminosity function of 
cluster galaxies. More important, we find that the relation L opt /Lx versus Lx is color dependent: Low X ray Luminosity clusters 
have a bluer galaxy population. 
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1. Introduction 

The estimate of the total cluster luminosity in different colors, 
the estimate of its Luminosity Function (LF) and the compar- 
ison with the LF of field galaxies (see also |Chincarini (19 88)) 
are powerful tools to gather information on how a cluster of 
galaxies forms, evolves and accretes from the field. 

Currently, we have at least two determinations of the form 
of the local field LF (Binggeli et al. 1988) which are in excel- 
lent agreement: [ Zucca et al. (1997)| and |Blanton et al. (200 1)| 
Blant on et al. (200 1)| estimate the local field LF in different 

colors down to an absolute magnitude M 16 (Petrosian 

magnitudes) using the Sloan Digital Sky Survey (SDSS) 
commissioning data. The very large sample they used al- 
lowed them to have very good statistics and they estimated: 
4> = (1.46 ± 0.12) ■ 10~ 2 h 3 Mpc- 3 , M*(r) = -20.83 ± 0.03, 
a — -1.20 + 0.03. These values are comparable with the 
ones found by |Zucca et al. (1997)1 wno analysed the local 

field LF down to magnitudes M 12 through the ESO 

Slice Project (ESP) survey: = (0.020 ± 0.004)h 3 Mpc~ 3 , 
M*(Bj) = -19.61 + 0.07, a = -1.22 ±0.06. The data from 
the ESP survey, four magnitudes deeper than SDSS, showed 
a steepening of the faint-end slope starting at magnitude 

Mbj 17 that Zucca et al. fitted with a power law with 

slope /3 = -1.6. In spite of the small number of objects with 
Mbj > -16, this is a very important result, as it is in line with 
the value found by |Driver^Uj^Q994)| and it is in contrast 
with the value that [LTveT?ay^Hu^^99"2")| obtained from a 



shallower survey, and highlights the need for further studies of 
the faint-end slope, at least down to magnitudes M 12. 

We believe that we have a fair knowledge of the local field 
LF, even though future surveys will provide better details, es- 
pecially at faint luminosities. Naturally, the number of low 
surface-brightness galaxies (LSB) we miss and the number of 
compact galaxies that could be confused with stars are still un- 
known. 

According to |Caldwell & Bothun (1987)] LSB galaxies are 
destroyed near the cluster center more likely than normal 
galaxies, and form a diffuse stellar background following 
the cluster potential. The compact galaxies, thanks to their 
strong density gradient, will survive also if their number is 
very high in a cluster, as found by jOrincari ni & Rood (1972)| 
IGavazzi et al. (2002)| and |Sakai et al. (2002)| Undoubtedly, the 
analysis of the complete SDSS will give an important step for- 
ward in solving these problems. 

For clusters of galaxies the situation is more com- 
plex. A milestone in the study of the LF is the cata- 
logue by Binggeli et al. (1985 a,b)| and the related analysis by 
|SandageetaL*(1985)| For the first time, they understood the 
role that the different types of galaxies play in the LF. They 
found that the faint end, clearly dominated by dwarf galaxies, 
had a rather high slope of a 1 .4. 

Because of the closeness of the Virgo cluster, these works 
should be taken as a firm reference point, and every study of the 
LF of cluster galaxies must be aware of the difficulties in de- 
tecting and measuring LSB galaxies (Impey & Bothun (1997)) . 
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A large amount of work has been published on the LF of 
cluster galaxies and we will discuss the recent literature in a 
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comparative way in a forthcoming paper. To give a short survey 
of the situation, we refer here to a few works. 

Studying A 963 (z ~ 0.2), |Driver etal. (1994)| find a re- 
sult that coincides very well with the early work on Coma 
(z ~ 0.02) by |Godwin & Peach (1977)] and with the LF func- 
tion derived for A 2554 (z ~ 0.01, |Smith et al. (1997)) . Their 
analysis on the sample, that reaches Mr = -16.5, shows a LF 
that could be fit by a composition of two Schechter functions, 
or by a Schechter plus a power law, with slope - 1 and -1.8 
respectively (in reasonable agreement with the findings from 
the Virgo cluster and by the ESP in the field). Also, this work 
outlines the uncertanties related to the background subtraction 
(Driver et al. (1994 Fig. 5a); see also |Bernstein et al. (1995)) . 
As also shown by Ab ell et al. (1989)| the background density 
of galaxies is uncertain by a factor of the order of the effect 
we measure and shows changes from region to region. From 
published studies and from our own work, this seems to be the 
main source of bias in determining the slope of the LF faint 
end. 

To better distinguish between the probable cluster mem- 
bers and the background, B iviano et al. (1995)| used the red- 
shift data in their work on Coma. Obviously, this would be 
the way to go, but generally redshifts are not available for 
the faintest objects. The spectroscopic sample used by Biviano 
et al., reaches magnitude Mb = -16.9 (assuming for Coma 
m - M = 34.9 with Ho = 75km/s/Mpc), with a completeness 
of about 95%. To extend the analysis to fainter magnitudes, 
M ~ 15.5, they added 205 galaxies selected by photometric cri- 
teria and the faint end slope they measure is about a 1.3. 

Further confirmation of those results was obtained by 
|Bernstein et al. (1995)| who used the deepest observations 
available of the Coma cluster (down to Mr = -9.4). They 
were not able to measure a detailed LF on the bright end, 
but found a slope a = -1.42 ± 0.05 in the magnitude range 
-19.4 < Mr < -11.4, and a sharp steepening of the LF be- 
yond that range. This probably recalls the old problem posed 
by Zwicky (1972 - private communication) about the faint end 
of the LF: where does it end and what kind of objects populate 
it. 

At faint magnitudes, it's indeed hard to distinguish the faint- 
magnitude objects of the cluster from background objects with- 
out high resolution imaging and spectra. 

Another interesting result has been pointed out by 
Bark house et al. (20 02) They analyzed a very large sample of 
cluster galaxies with the same purpose as ours, that is to de- 
termine if the faint-end slope is a function of the cluster mor- 
phology and if there is a gradient in the faint end slope of the 
LF moving from the central cluster regions toward the out- 
skirt. They seem to have reached some positive evidence on 
those effects even if they do not go very deep in magnitudes 
(Mr < -16). Finally, analysing the LF of the first cluster of our 
sample, A496, Molinari et al. (1998), found a clear evidence of 
bimodal LF and a quite high faint end slope (a = -1.65). 



2. The luminosity Function 

The detailed description of the observations, the data analy- 
sis and the algorithms we used can be found in the papers by 



|Molinari et al. (1998)||Moretti et al. (1999)| and the laur ea the- 
sis by |Moretti (1997)||Ratti (1998)| and |Parolin (2002)| As we 
want to present the results from the analysis of 7 clusters of the 
original sample, we give here just few informations to ease the 
reading (see TableQ. 

We carried out the observations at ESO, La Silla, at 
the Danish 1.54 m Telescope using the Danish Faint Object 
Spectrograph and Camera (DFOSC) and the g, r and i col- 
ors of the Gunn's photometric system (T huan & Gunn 19761 
IWade et al. 19791 . 

After a standard reduction of the data, for each cluster 
we created a catalogue of sources by merging the photo- 
metrical catalogue of the bright and extensive galaxies ob- 
tained through a program specifically developed by Moretti et 
al.(1999) with the photometrical catalogue of small and faint 
objects obtained using the MIDAS Inventory package. Infact, 
as Inventory was developed to detect sources in distant clus- 
ters (point-like sources), we first extracted and analysed all the 
bright and extended galaxies through the program by Moretti 
et al. (1999), that reproduces the image of a galaxy through the 
Fourier analysis of its isophotes and subtracts it to the main 
frame, then applied Inventory. 

The catalogs thus created were tested through a series of sim- 
ulations, briefly described below, to evaluate their complete- 
ness: using one of the galaxies reproduced with the program 
by Moretti et al. divided by coefficients ad hoc, we created a 
sample of artificial galaxies with a known distribution of mag- 
nitudes. We partitioned the field to be tested in circular annuli 
in which we added the artificial sample with known random 
coordinates; by running Inventory on the added fields we were 
able to estimate the selection function as a function of the dis- 
tance of the cluster center. 

After that, we carried out our analysis. 

2.1. The standard analysis 

To sharpen the contrast between the counts in the cluster and 
the counts in the field, we selected those galaxies that in the 
color magnitude plot are placed in between the 68% confidence 
curves defined by the line fitting the r versus (g-r) distribu- 
tion of bright galaxies (the equivalent of the E/S0 sequence as 
defined by |Visvanatan & Sandage (1977)| l and its rms (Fig.Q. 
These galaxies, mostly ellipticals, have indeed a higher proba- 
bility to be cluster members. 

This area is defined by the curves: 



C = a- m + q±(k- crfl t + ctq) 



(1) 



Where the term a • m + q represents the fit of the color- 
magnitude relation for galaxies brighter than r = 19 and 
o"o = VT^g + of) is the color index error derived using the 
magnitude error estimated on the fields overlapping in the dif- 
ferent frames. 

Non-members may perturb the fit to derive the coefficients a 
and q, also in view of the small number of stars which belongs 
to the bright sequence, so we coupled the regression solution 
with the best fit as estimated by eye. 

We assumed that all the galaxies which in the color- 
magnitude diagram belong to this region form what we will 



I, Parolin et al.: The Luminosity Function of Cluster Galaxies. III. 



3 



Cluster 


Q"2000 


(52000 


Redshift 


B.M. type 


Abell type 


Abell 0085 


00h41m50.11s 


-09dl8ml7.5s 


0.05560 


I 


2 


Abell 0133 


01h02m42.21s 


-21d52m43.5s 


0.05660 


I 


2 


EXO 0422-086 


04h25m5 1.02s 


-08d33m38.5s 


0.03971 


III 


- 


Abell 3667 


20hl2m35.08s 


-56d50m30.5s 


0.05560 


II 


2 


Abell 3695 


20h34m46.86s 


-35d49m07.5s 


0.08930 


I 


2 


Abell 4038 


23h47m41.78s 


-28d08m26.5s 


0.02920 


III 


2 


Abell 4059 


23h57m00.02s 


-34d45m24.5s 


0.04600 


I 


1 



Table 1. The 7 clusters analized in this paper. 
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Fig. 1. A0085. The total counts over the whole area observed 
are plotted together with the background counts as measured 
by |Tyson (1992)| Due to the large background component on 
the outermost area for the faintest galaxies, the counts on the 
cluster galaxies almost coincide with the background counts. 



define sequence galaxies and charachterize the cluster sample. 
These are the objects we used to estimate the parameters of 
the distribution function by the method of the maximum likeli- 
hood. 



Cluster 


All 


Sequence 


Sequence 


Rnax 








(within 400") 


(arcsec) 


A0085 


2701 


1688 


659 


1473.2 


A0133 


4022 


3068 


1016 


1350.6 


A3667 


3205 


1887 


521 


1597.9 


A3695 


2156 


1429 


751 


1293.1 


A4038 


3955 


2949 


767 


1714.3 


A4059 


1778 


1201 


543 


1107.0 


EXO 0422-086 


1612 


1254 


604 


1157.8 



Table 2. Number of galaxies for each cluster and maximum 
distance from the cluster center covered by the data (R max )- The 
clusters have been observed up to an everage distance of about 
23 arcmin or 1 .08 Mpc, that is the observations cover almost an 
Abell radius and therefore we observed up to a distance from 
the center of the order of r2oo- 



We assumed that the sequence galaxies are, as stated above, 
in large part cluster galaxies. This is true especially for the 
brighter objects where the background is pratically absent. In 
other words, we subtracted only those objects outside the se- 
quence where it is more likely that the background dominates 
at faint magnitudes. 

While this assumptions has no effect in the estimation of the 
Luminosity of the clusters since that is due essentially only to 
the bright cluster objects, it may seriously affect the estimate 
of the slope of the faint end of the LF. Indeed instead of mea- 
suring the slope of the faint end of the Luminosity Function 
we could be measuring the slope of the background counts. To 
demonstrate that by selecting the sequence galaxies as the sam- 
ple on which to measure the LF parameters we do not bias the 
background contamination, we estimate, for a few particularly 
rich clusters, the Luminosity function after background sub- 
traction and obtained values in perfect agreement with what 
we obtained by the analysis of the sequence galaxies and the 
use of the maximum likelihood algorithm. 

The background counts of galaxies in the Gunn and Thuan 
photometric system were estimated using the counts published 
by |Tyson (1992)] 

The background counts as measured by |Tyson (1992)| were 
used instead of the SDSS counts since, at the time we started 
the observations and the analysis these were, at the best of our 
knowledge, the best counts we could use at different wave- 
lengths. Furthermore, these counts were in good agreement 
with the counts measured in our fields using the outermost 
fields observed in those clusters where we estimated the back- 
ground. Moreover the SDSS counts cover a smaller dynamical 
range and stops at m ~ 20 and the early attempts we made, 
using preliminary SDSS data as distributed in the Web, to ex- 
trapolate to fainter magnitude the slope (down to about r ~ 25) 
showed that we were overestimating the background counts. 

Indeed by plotting the coefficients of the logarithmic fits of 
the counts given by Tyson as a function of the effective wave- 
length of the pass-band filter, we notice that these observed val- 
ues are accurately fitted by a polynomial function. It is then 
straightforward to estimate by interpolation the parameters ex- 
pected for counts in the Gunn and Thuan photometric system, 
by so doing we obtain the following relations: 



logN g = (0.45: 
log N r = (0.39 : 

logN; = (0.34: 



: 0.056)g - (6.30 ± 0.89) 
0.063)r - (4.67 ± 1.01) 
0.070)i - (3.38 ± 1.11) 



(2) 



The errors given in the above relations represent the pro- 
gressive errors estimated on the interpolation relation by prop- 
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agating the errors of the fits relating the parameters to the ef- 
fective wavelength of the filters used by Tyson. 



Fig. 2. A0085. Counts in the r filter limited to the central frame 
compared to the background counts as given by |Tyson (1992)| 
The cluster signal is well visible over the background. 

We modeled the distribution in luminosity by using a bi- 
variate LF that is the sum of a Gaussian and a Schechter. 
The choice is the most reasonable one following the work by 
ISandage et al. (1985)| on the Virgo cluster; it is evident from the 
data of the present analysis, Fig.^and Fig|2]and as indicated 
also in Molinari et al. (1998)| and references therein. In a clus- 
ter, as it is very evident also from our data, the bright galaxies 
are located in the central region and these are the galaxies that 
are fitted by the Gaussian component of the LF. That is the two 
component fit is valid only for the central region, in the out- 
skirts of the cluster the luminosity of the galaxies is distributed 
according only to the Schechter function (Fig. [8}- 

We are aware that not all of the published analysis on the LF 
of clusters of galaxies, such as the one of G oto et al. (2002)| ev- 
idence the bimodal distribution. While it is hard to answer why 
without remaking the analysis and eventually the data reduc- 
tion, it seems to us that the main reasons reside on the statistics, 
accuracy of the magnitudes and background subtraction and, in 
part, from the binning. Indeed from a practical point of view, 
and looking at it in a slightly different way, the main evidence 
depends from the small "gaps" in the observed distribution 
which is evidenced in the point where the composite functions 
are both somewhat weaker. In our cluster sample this occurs at 
the r mag which is in the range 17-18. Obviously in a redshift 
sample the signal vs. noise ratio is large due to a better selec- 
tion of cluster members detaching the cluster galaxies almost 
completely from the background. Indeed this is clearly shown 
also by the Coma sample selected by | Biviano et al. (1995)| and 
by the sample in More tti et al. (1999)| On the other hand the 
analysis o f|Goto et al. (20 02) based on the SDSS survey clearly 
shows that a) the analysis is consistent with the distribution of 
two underlying populations, b) the analysis deals with the com- 
posite LF so that details are smoothed out adding the contribu- 



tion of different clusters and c) the limiting magnitude is not 

faint enough (the sample stops at M 18) to clearly define 

the faint end. 

For the Maximum Likelihood fitting we used only objects 
brighter than the 14th magnitude. This means that we excluded 
the cD (in two clusters this cut off excludes 2 or 3 galaxies) 
from the fitting. The cD does not necessarily conform to a dis- 
tribution function as the other galaxies and its luminosity may 
be a function not only of the cluster mass distribution at the 
time of formation but, above all, may reflect the characteris- 
tics of evolution. While it is easy to account for its luminosity, 
it is a perturbing factor in the fitting. The cD galaxies will be 
discussed in detail in a forthcoming paper. 

We analyzed all the clusters in a similar way to have inter- 
nal consistency. In most cases, two of the authors analyzed the 
clusters and compared the results to have an estimate of how 
much the result could depend from our method of analysis and 
have, at the same time, an estimate for the uncertainties. The re- 
sults of the standard analysis are listed in Tables[3]and|3] Using 
these values and the master photometric catalogues we derived 
the Luminosities listed in Table[5] 



3. Checking our analysis, a few test cases 

In this section, we illustrate part of the data and our analysis 
also evidencing some weak points. While the standard analy- 
sis was done in a very homogeneous and systematic way, for 
most of the clusters we used alternative analysis, that is a) we 
estimated the luminosity function subtracting from the central 
frame the counts as given by |Tyson (1992)| b) we subtracted 
the background using the outermost field we observed where 
the cluster galaxies contamination is known to be negligible. 
In all these cases, after the subtraction of the background, the 
LF parameters were estimated after binning the derived his- 
togram, assumed to represent statistically the counts of the 
cluster galaxies, and using a% 2 minimizaton program. 

As it is explained below, it makes a difference whether or 
not the cD galaxy is part of the sample. Furthermore, and as 
expected, if the analysis is not limited to the central region, of 
the order of a core radius, the contrast cluster vs. background 
is largely washed out. But, as described in the following text, 
the best evidence that the bimodal distribution is a real fact is 
illustrated in the lower panel of Fig. [5] Here the contamination 
of the background galaxies brighter than the 21 st magnitude 
is negligible and the maximum of the Gaussian distribution is 
readily visible at r ~ 16 while it is also clearly visible the bright 
end of the LF and to the low number of galaxies expected to- 
ward the faint end of the Gaussian distribution. 

Since all the set of clusters used to test the method were 
treated in a similar way, for each cluster we discuss only part 
of it in order to give a complete view of the tests we made in 
oder to develop a feeling for the robustness of our results. By 
doing so we also avoid useless duplications. We also definitely 
think that the standard analysis and the results we give in Tables 
[3] and |4] represent the consistent results, with their errors, that 
are the output of this work and that is why we prefer not to 
confuse the issue by listing the derivation of the parameters by 
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other methods that while in agreement within errors with the 
standard analysis are less homogeneous and rigorous. 

3.1. A0085 

In Fig.^ we plot the total counts, that is the counts made over 
the whole observed area of our sample, for the 3 colors used 
(g, r and i), as a function of the magnitudes. Moreover, we 
plot the counts for the field as derived using the analysis by 
|Tyson (19 92) in the filter r and i. 

Our counts do not differ much from the background from one 
filter to the other. The larger discrepancy is with the i filter in 
which we count more than 30% of galaxies less than in the 
other filters. A probable cause for this is the high sky back- 
ground and the lower quality of these frames. However, we 
want to evidence the excess of counts (bump) compared to the 
background at magnitude ~ 15 that is due to the cluster we ob- 
served. 

If we plot only the central field, the contrast cluster/background 
incerases considerably and it evidences indeed the fact that the 
cluster galaxy counts can be clearly separated. This is naturally 
what we expect. Assuming a King's density profile the percent- 
age of cluster galaxies we would expect on the outermost field 
we observed in the different clusters is of about 10% so that 
here the backgrounds dominates. Indeed, and always to check 
out procedures, we also estimated the cluster galaxies counts 
by subtracting the outermost field from the central once assum- 
ing that the small cluster contamination on the background so 
defined would not bias the estimate od the cluster LF. We got 
a good agreement, for those clusters which were analysed also 
in this way, with the standard analysis. 

The bright end is most highly dominated by cluster galaxies 
as we expect a small number of field galaxies in such a limited 
volume of space. 

After the correction for completeness, it is straightforward 
to estimate the cluster contribution by subtracting the field 
counts. However, by doing so, we are more sensitive to the 
fluctuations of the background and to an eventual uncontrolled 
incompleteness introduced by our analysis. Furthermore, at the 
bright end we are dealing with a few objects, and the low statis- 
tics make things more uncertain. That is why we decided to use 
only relative measurements and, to further increase the contrast 
cluster versus field, to use only a selected subsample of galaxies 
for each cluster delimited by the fitting of the bright sequence 
galaxies and the lcr error. 

As we explained earlier, we define a color-magnitude se- 
quence via the color magnitude diagram for each filter and for 
each cluster. In the case of the r observations for the cluster 
A0085, Fig.[0 the lines delimiting the sequence galaxies area 
are defined by: 

(g-r) S eq = -0.0180-r + 0.802 

crerr-cur = (g - r) seq + (0.1 17 + 0.675 ■ 0.217) 

For this cluster, we used a factor 0.675 to reduce to 50% the 
probability that a galaxy belongs to the sequence with the com- 
puted rms. We also extended the sequence multiplying by a fac- 
tor 3 the value of cr^t to increase the number of sources at the 
faint end and to ease the comparison with the other clusters (in 




-1.5 -1 -0.5 0.5 1 1.5 2 



Fig. 3. The color magnitude plot for the central area of the clus- 
ter A0085. The plot illustrates the fit, shown here as a dashed 
black line, that determines the sequence galaxies as described 
in the text. The solid lines, as described in the text, define 68% 
probability that the objects has a color-magnitude as defined by 
the sequence. 

particular A4059). These minor adjustments do not affect the 
derivations of the parameter in which we are interested. On the 
other hand, they are selected during the analysis to optimize the 
statistics and the homogeneity of the solution. Two members 
of our group derived the luminosity function both using the se- 
quence we just defined and all of the central field objects. This 
checking analysis was also carried out using the;f 2 method and 
subtracting the background. For the faint end of the LF, the re- 
sulting slopes were all very similar to that obtained with the 
standard analysis. In all cases, we obtained similar and consis- 
tent results within the errors, an evidence that the results are 
good and practically unaffected by the method of the analysis. 

3.2. A0133 

With this cluster we had some difficulties in getting a reason- 
able fit of the observed LF in spite of the large number of mea- 
surements. As we did for most of the clusters, we carried out 
the fits in the different wavelengths deriving the parameters of 
the LF: 

a) for all the galaxies in the central field after subtracting the 
background; 

b) for the sequence galaxies of the cluster after background 
subtraction; 

c) for the counts derived as the difference between the central 
field and the outermost field (for the clusters for which we 
could use the latter as a background). 

In spite of the large uncertainties, we derived the same parame- 
ters within errors with all methods. Finally, the various samples 
were tested for solution in a semi-empirical way. The counts 
as a function of the magnitudes, were fitted by modifying the 
parameters derived statistically and estimating the best fit, or 
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Fig. 4. Color-color plot of A0085. The bright sequence galaxies 
lays exactly on the color-color curve generated by the models 
of synthesis. Naturally, the sub-sample sequence fits as well 
with a broader dispersion and reflects the color-color cuts we 
imposed in the definition of the sequence. 



the different fits, also by eye. This procedure, that we describe 
only for A0133, was indeed carried out for most of the sample 
clusters. The reason is that in some cases the parameters we 
derived formally using the Maximum Likelihood have large er- 
rors. This is especially true for the magnitude of the knee of the 
Schechter function or for the ratio of the Gauss to the Schechter 
normalization factor. By doing so we hoped to avoid flukes due 
to the analysis and to the small counts and to make sure which 
parameters must be taken with caution. What remains certain, 
however, is that the LF is in all cases dominated by a rather 
broad Gaussian defined by the brightest galaxies and dominat- 
ing over the Schechter function which is essentially defined by 
the fainter galaxies. The cD galaxy does not fit the Luminosity 
Function. 
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Fig. 5. Abell 4038. Upper panel: the completeness Function for 
the r filter as a function of the distance from the cluster cen- 
ter (in seconds of arc). Lower panel: the distribution of the se- 
quence galaxies shows that the number density, and therefore 
the distribution function we derived, is composed by a con- 
densation of bright objects and by a multitude of faints objects 
whose number increases with the magnitude. 



3.3. A4038 

With a completeness down to the 22nd magnitude in the central 
region and down to the 23rd magnitude in the outskirts (upper 
panel of Fig.|5}, A4038 is a clear example of the bimodal dis- 
tribution function. 

The lower panel of Fig. [5] shows the color-magnitude plot 
of the sequence galaxies. We see quite clearly the increase first 
and the decrease later of the number density of galaxies along 
the sequence going from brighter to fainter magnitudes, and 
this effect has been observed in all clusters and in all the col- 
ors. At about r = 19, we find the so called gap, that is the region 
in the LF where the Gauss distribution dominating at bright 
magnitudes merges with, and is taken over, by the Schechter 
distribution. 

In the outermost field of this cluster the background, as esti- 
mated by Tys on (1992)| is somewhat higher than our counts for 
r > 21 .5. This outlines that at faint magnitudes we are strongly 



affected by the fluctuations of the background and by incom- 
pleteness. 

3.4. A4059 

We use of this cluster to discuss in more detail the fit and related 
uncertainties on the study of the LF of clusters. These uncer- 
tainties are not due to pure statistics, but rather to the selected 
procedures. As usual, we illustrate the r filter observations for 
uniformity with the previous, and following,discussion. A sim- 
ilar analysis has been done and checked in all colors on vari- 
ous clusters obtaining consistent results. That is to say that the 
goodness of the parameters values obtained and the procedures 
discussed is robust. 

In Fig.|6j we plot both the data obtained by subtracting the 
background as estimated by |Tyson (1992)| and that derived by 
our data using the outermost field observed (as we stated earlier 
here we have a cluster galaxies contamination of about 50%). 
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Fig. 6. A4059. The black filled diamonds refer to the counts 
in the central 400" frame after subtracting the counts by 
|Tyson (19 92) while the empty dots refer to the same counts 
after subtraction of the outermost field we observed for this 
cluster. The solid line reflects the values derived from a semi 
empirical fit of the data. The dash line show the Gaussian dis- 
tribution selected and the fit of the faint end by a straight (dot) 
line. The fit has been obtained using the following values for 
the LF parameters: G/S = 350.0/150.0, m G = 16.3, cr c = 1.0, 
ms = 18.0 and a = -1.5. The cD galaxy does not fit the LF. 



The Tyson Background is somewhat higher, the difference be- 
ing probably due to fluctuations. Contrary to what has been 
done for standard solution, in the distribution in Fig.[6]we also 
included the cD galaxy that, being too bright, does not fit the 
LF, no matter which parameters we select. 

At the faint end we obtain a fit of the counts that is in ex- 
cellent agreement with the value obtained using the Maximum 
Likelihood method applied to the sequence galaxies. This jus- 
tifies the method once more and underlines that we derive re- 
liable values also for the faint end in spite of correcting only 
partially for the background galaxies. 

However, to stress the subtleness of the faint end fitting 
in presence of a background contamination, we point out that 
this slope can easily be estimated by a linear fit of the faint 
galaxies. The faint end of the Schechter function expressed 
in magnitude and on a logarithmic scale is given by the re- 
lation -0.4 ■ (m - m*) • (or + 1) so that a can be derived by a 
simple linear fit of the faint end of the observed distribution 
using the relation: a = - (1 + (observed slope /0.4)). The back- 
ground can also readily checked and subtracted since we sim- 
ply make the difference between two straight lines. Assuming 
the counts are dominated by background objects, we would 
have observed slope = 0.39 and derive a > -1.5 and closer to 
-2. 



4. The cluster luminosity: the contributions of the 
bright and the faint end 

An estimate of the luminosity of a cluster is given by simply 
integrating the Cluster LF, in absolute magnitudes, using the 
following expressions: 



L s = 
L G = 



r % / , \a+l 



1.086 



plogL 2 
JlogLl 



-0.5 



-2.5(loeL-lQ8L G ) ¥ 

M(logL) 



(3) 



where Li > L2 and where rmp is the Gaussian/Schechter nor- 
malization ratio. 

As we said before, in the expression of the LF the faint end 
is the greater source of uncertainties because of the background 
fluctuations and of the uncertainties in the estimate of the com- 
pleteness (the estimate of the selection function remains often 
uncertain in spite of the Monte Carlo approach which is less ro- 
bust when applied near the core area where many bright galax- 
ies are located). 




Fig. 7. The Ratio between the luminosity contained in bright 
galaxies fitted by a Gaussian LF and the total luminosity 
(Gaussian + Schechter function) as a function of the param- 
eter a and the normalization ratio between the Gaussian and 
the Schechter function, p in the equation. From top to bottom 
p = 2.0, 1.0, 0.7, 0.4, 0.2. 

It seems that the cluster luminosity is dominated by the 
bright galaxies. In Fig. 0for various values of the ratio be- 
tween the Gaussian and the Schechter normalization, we plot 
the percentage of the total cluster luminosity which is due to 
the bright galaxies as fitted by a Gaussian. We used a magni- 
tude difference < Mcauss > - MJj chechtei . = 1.5. This difference is 
close to what we find in some of our clusters. In the Virgo clus- 
ter, |Sandage et al. (1985)1 find < M Gauss > -M; chechter = 0.9. 
However, the ratio is practically not affected very much by such 
differences. As it can be seen from Fig.0 what counts is mainly 
the normalization ratio and the value of the faint end slope a. 
Disregarding extreme values, that is p > 1, which have been 
however also observed in our clusters, we find that for a 1 .4 
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the ratio varies between 0.6 and 0.85 while for a 1 the 

ratio changes between 0.74 and 0.92. Clearly, in those cases 
where p > 1, the Gaussian distribution function contribution 
dominates. In the Virgo cluster (Sandage et al. 1985) we esti- 
mate p ~ 0.14. However, in this cluster <x G is rather large while 

a 1.35 close to what has been found in the majority of the 

clusters. 

Merging clusters, or more precisely clusters for which sub- 
structures are detected, do not affect the estimate of the total 
luminosity and the LF. During the merging of substructures we 
expect changes in luminosity of the galaxies due to induced 
tides and therefore luminosities for the merging galaxies that 
differ somewhat for the luminosities of the merged galaxies. 
On the other hand this effect, which likely result in a loss of 
stars to the Intra-Cluster Medium (ICM) and at the same time 
brightening and star formation due to the induced turbolence in 
the Interstellar Medium (ISM), is a secondary, not detectable 
effect to the analysis of this work. 

We calculated that an estimate of the cluster luminosity us- 
ing only the bright galaxies as fitted by a Gaussian plus a stan- 
dard slope faint end luminosity function gives a mean error of 
about 10%. In all those cases in which there are uncertainties 
in the faint end fitting or we notice the presence of a greater 
background contamination, we estimated, as a check, also the 
total luminosity by using mean parameters for the faint end of 
theLF. 

In general, the procedure we followed is the following: 

- We used the cluster parameters as derived from the standard 
analysis 

- We normalized the luminosity function by subtracting the 
luminosity derived by the counts of the background as es- 
timated by |Tyson (1992)| from the counts we have for each 
cluster in a given range of magnitudes selected ad hoc (the 
range is not critical): g < 19.5, r < 19 and i < 19. 

- We limited the normalization to the bright galaxies in order 
to avoid further uncertainties due to the background and the 
subtraction of large numbers at faint magnitudes. 

- In most cases, we based the normalization on the Gaussian 
component since the ratio of the Gaussian to Schechter was 
estimated from the fit, albeit with large error in some cases. 
For the normalization, we used only the frame about the 
cluster center. 

- For each cluster we also counted and selected the galaxies 
in the range of apparent magnitude corresponding to M g < 
-16.5, M r < -17.0 and M,- < -17.0 and we simply added 
up the relative luminosities. 

We observe a clear difference between the central part of 
the cluster, where the bright galaxies dominate, and the out- 
skirts, where the central part of the cluster bright galaxies are 
generally absent. This is a strong evidence of segregation in lu- 
minosity and is consistent with the model of a fairly relaxed 
cluster of galaxies. That has been taken into account in com- 
puting the cluster luminosity. 

To better account for this difference we adopted a model in 
which the number density distribution follows a King's profile 
and the LF changes from a Gaussian plus a Schechter (where 



the luminosity of the Gaussian is dominant) to 
function from the central regions to the outskirt the 



133 - early-type R < 400 



a Schechter 
cluster (Fig. 




Abell 133 - early-type R > 1200 




18 20 
mag r 



Fig. 8. A0133. For the sequence galaxies within 400" from the 
cluster center the LF is a Gaussian plus a Schechter (upper 
panel), while for the sequence galaxies more distant than 1200" 
from the cluster center (lower panel) the LF is well fitted by a 
Schechter function. 

To estimate the contribution of the outskirt of the cluster to 
the luminosity, we thus assume a decreasing density obeying 
the King's law and account for the fact that we do not observe 
bright galaxies in the most distant regions of the cluster. We ex- 
trapolate the estimate to 1 Abell radius (AR = 1.5h _1 Mpc) and 
we use a core radius of 0.25IT 1 Mpc in all those cases for which 
we do not have a direct estimate of the core. The normalization 
of the King's function, 0*(r), is given by the relation: 



£ r(r)0(M)drdM = JT 



tl.Field 



> 0(M)dM 



(4) 



with r e 



Area(Cl.Field) 



where with "CI. Field" we generally 



refer, unless otherwise stated, to the observed central frame. 

As we have shown in the previous section, the cD galaxy, 
or the two brightest galaxies, hardly fit the analytical functions 
we used. Indeed depending on how the brightest galaxies form, 
or grow, it may be reasonable to expect a peculiar behavior and 
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disagreement with the analytical function. In other words the 
analytical fit may underestimate the Luminosity because such 
a bright galaxy is not expected by the distribution. To com- 
pute the luminosity, the LF has been considered in the range 
between M = -12 and M = -24. 

4.0.1. A0085 

To have reasonable statistics and a small background contami- 
nation, we normalized the luminosity function with the counts 
of the cluster galaxies within a distance of 300" from the cluster 
center, that is about the area of the cluster core. The observed 
area was computed from the geometry of the observations and 
accounting in an approximate way for the occultation due to 
the cD and a few other bright galaxies (a rather small correc- 
tion however). To avoid a large correction for incompleteness, 
we also limited the counts to galaxies brighter than r = 23.0 
(M r = -14.61 + 5 log h), while for the background we used the 
counts as given by T yson (1992)| 

Normalizing, the luminosity distribution function using the 
observed galaxies within 300", we compute a luminosity in the 
r band of L 30 o<< = 9.1 • 1O 1O L , with 46% of this given by the 
bright Gaussian component. The contribution of the Schechter 
component within one AR is Lar = 1.08 • 1O U L , so that the 
total luminosity of the cluster (as we said there are no bright 
galaxies outside the core radius) within an Abell Radius is 
LibtAR = 1-99 ■ 10 n L G . A somewhat larger value is obtained if 
we normalize the LF to the bright galaxies only. In this case, 
we derive L300'' = 1.26 • 1O U L and L To1 ar = 2.76 • 1O U L . 

The parameters we derived using the standard solution for 
the LF, Tables [3] and @] are rather stable in the sense that were 
obtained, in agreement within errors, using different methods: 
X 1 and Maximum likelihood. Unlike most of the other clus- 
ters we analysed, we find a value of p (the ratio between the 
Gaussian and the Schechter function coefficient) larger than 
one, a rather small era and a magnitude for the knee of the 
Schechter function (m*) which is brighter than the mean mag- 
nitude of the Gaussian distribution (mc). This is suspicious and 
certainly due to the fluctuations related to the very few bright- 
est galaxies. We can account of this by using a semi automatic 
(in part also by eye) fitting with the result, quite reasonable, of 
the lower panel of Fig.|5J where the brightest galaxy has also 
been plotted. 

4.0.2. A0133 

For A0133 we measure a rather steep faint slope. However, 
we may be somewhat contaminated by background on 
the faint end. By normalizing on the bright galaxies 
(12<r<19) detected in the central frame (this has a ra- 
dius of about 0.27IT 1 Mpc) we derive a Gaussian luminos- 
ity Lg = 1.73 • 10 n L o and a total luminosity within an Abell 
Radius LjotAR = 2.16-1O U L . Due to the small contribution 
by the Schechter Luminosity (for all clusters the Luminosity 
is given in the magnitude range M = -24 to M = -12) by us- 
ing a = -1.4 rather than a = - 1 .65 the luminosity decreases of 
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r mag 
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12 13 14 15 16 17 18 19 20 21 22 23 

Fig. 9. Two possible fits of the LF for A0085: excluding the cD 
galaxy (upper panel), where the errors are simply the square 
root of the counts, and with a different binning trying to include 
with an optimization by eye the cD galaxy (lower panel). See 
however the appendix for the formal solution without cD. 

7%. A smaller value, however, is derived using the maximum 
likelihood fit. We adopt L Tot AR = 2.16 • 10 u L o . 

4.0.3. A4038 

The luminosity of A4038 has been computed using the 
fit to the bright objects, while for the parameter of the 
Schechter function component we used the mean, when ap- 
plicable, of the values estimated for the three colors by us- 
ing the Maximum Likelihood function, that is: a = -1.39, 
M*(r) = -13.81 + 5 log(h). For the ratio between the Gaussian 
and Schechter part of the LF, we compute G/S — 0.14. This 
value could be strongly dependent on the contamination due to 
background galaxies, so that we give as a reference also the lu- 
minosity computed for G/S = 0.8 in order to estimate the max- 
imum uncertainties we may have in the result. For the Gaussian 
fit we use Norm = 300, mc = 16.0 and <x G = 1.0. 

In the range 12 < r < 19, that is -22.7 < M r < -15.7, 
we observe 30 galaxies in the central field of the cluster. 
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r mag 

Fig. 10. A0133. Black dots show that the separation between 
the Gaussian and the Schechter is not that clear and even- 
tually could also be fitted, with similar uncertainties, with a 
Schechter. 

Subtracting the galaxies in the outskirt of the cluster we have 
a total of 26, while subtracting the background counts as esti- 
mate by Tyson (1992) we are left with 20 galaxies. These dif- 
ferences are minors and affect the computation of the cluster 
luminosity in the central field of only a few percent, so that 
they could be easily disregarded with respect to other uncer- 
tainties. For normalization, we use Nobs = 26. We also tested 
for uncertainties due to the estimate of the magnitude of M* 
in the Schechter LF. Field contamination would indeed make 
this estimate uncertain and favor fainter magnitudes. The to- 
tal luminosity changes by at most 5% by brightening M* of 
about 2 magnitudes, these value are given below in parenthe- 
sis. We conclude that the estimate of the cluster luminosity, 
dominated by the bright galaxies, is rather robust. The estimate 
gives: L c i ustel - = 6.7 • 10 10 (7.8 • 10 10 ) L with a contribution by 
the bright galaxies of 98% (84%). If we change the ratio the two 
Luminosity Functions, that is making G/S = 0.5 for instance, 
the contribution of the Schechter part is reduced to about 5% 
while the contribution of the bright galaxies remains constant. 

For this cluster we normalized on the number of objects 
that have a distance from the cluster center that is smaller 
than 400". Using the adopted parameters, we estimate that the 
luminosity due to faint galaxies within an Hubble radius is 
Lfaint = 3.7 • 1O 9 L . Since there is no additional contribution 
by bright galaxies, the toatal luminosity of the cluster would 
be Ljotai = 7.0 • 10 10 L o . Assuming the knee of the Schechter 
LF to be 2 magnitudes brighter, we derive a somewhat brighter 
luminosity, Lj tai = 9.8 • 10 10 L o . In this case, however, the 
number of expected galaxies would have been larger (525 
rather than 139) and we would have detected this larger over 
density. We adopt L Tot AR = 7.0 ■ 10 I() L o . 
The cluster shows three galaxies brighter than apparent 
magnitude 14.0 which do not differ much in luminosity. 
The brightest galaxy has a luminosity L g = 0.67 ■ 10 n L o , 
L r = 0.51 • 1O U L , L; = 0.39 • 10 n L o . By adding the 
three brightest galaxies we have: L g = 2.12 • 10 n L o , 




13 14 15 16 17 18 

r mag 



Fig. 11. A4038. The fit using the counts (black dots) obtained 
subtracting the outskirt cluster region from the central region 
gives a normalization value of 350 counts per magnitude per 
square degree. The counts obtained by subtracting the back- 
ground as estimated by |Tyson (1992)| diamonds, give a some- 
what smaller value. We adopt a normalization factor equal to 
300. 

L r = 1.68 ■ 1O U L , Li = 1.56 • 10 u L o . And this clearly 
shows: 

a) that most of the light is contained in the central bright 
galaxies; 

b) that sometimes we can not measure the cluster luminosity 
simply by fitting the counts. 

4.0.4. A4059 

As discussed in Sec. 13.41 the different estimates of the param- 
eters of the LF for A4059 are not in perfect agreement and de- 
pend somewhat of the method of analysis. We adopted those, 
as we illustrated there, which better satisfy the visual inspec- 
tion and, in spite of smaller statistics, have been derived after 
careful sky subtraction. Also in this case the cD does not to 
fit properly the luminosity distribution function: its luminosity 
indeed is larger than the whole Gaussian contribution so that 
we should correct for this later on. The luminosity of the cD 
is L g = 1.72 ■ 1O U L , L r = 1.88 • 1O U L , L, = 1.56 • 1O U L . 
The total cluster luminosity is Ljoiar ~ 1.0 ■ 1O U L/L and 
LTotAR = 1.9 ■ 1O H L/L applying a correction for the cD. 

The luminosity of the core without accounting for the 
cD galaxy and using both the adopted distribution func- 
tion and the others as derived in the respective the- 
sis work are: L core = 7.5 • 1O 1O L ; L^ iK - = 4.9 • 1O 1O L ; 



4.1. The luminosity within 400" from the cluster center 

To check for consistency with the total luminosity computed 
using the parameters of the LF and the King's profile, we also 
measured the luminosity of the central field of each cluster, 
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that is within 400" from the cluster center of the observed 
area, due to galaxies brighter than M g = -16.5, M r = -17.0, 
and M; = -17.0 by simply adding up the luminosity of the sin- 
gle galaxies. To the luminosity so computed, we subtracted the 
contribution of the background, Eq.|2 with: 



N„ =5128 , N r = 6680 , Nj = 7920 



/ 

*j in 



and the luminosity of the background computed as: 

.m(M„=-16.5,M„i=-17) 
m=11.5 



N(m)L(m)dm 



(5) 



(6) 



since we did not detect cluster galaxies brighter that m = 
1 1 .5 in any filter. 

Only for the cluster EXO 0422-086 we used magnitude cor- 
rected for galactic extinction. The cluster magnitudes estimated 
in this way are not affected by any bias due to the analysis. 
However, more distant clusters (as for instance A3695) will 
appear brighter since a larger cluster area has been accounted 
for. The effect is however rather small, and we can account for 
it anyway, for we are considering only fairly bright galaxies. 
Trying to have results as robust as possible, we computed also 
the Cluster Luminosity as follows: using the estimate of the LF 
derived by the Maximum Likelihood method we computed the 
luminosity of the central region observed, Field 1, normaliz- 
ing to the number of galaxies observed in this Field. The Field 
covers the observed galaxies within R < 400 arcs and, except 
for the cluster A3695, is of the order of the core radius. We 
add the contribution of the brightest galaxies excluded from 
the fit: only one with m < 14.0 for all the clusters but A3695, 
for which the fit was carried out for galaxies m > 15.5. To com- 
pute the cluster luminosity to the Abell radius (1.5/i~' Mpc), as 
we indicated before, we add the contribution of the Schechter 
function component weighted by a King's profile for account- 
ing of the density distribution. We used a core radius of about 
0.25 Mpc except in the case of A0085 and A4059, clusters for 
which we were able to measure a core radius, albeit with rather 
large errors. For each cluster, we subtracted from the galaxy 
counts the background interpolating for the filters we used the 
counts given by |Tyson (19 92) in different colors. 

As it was mentioned before, the Maximun likelihood fit 
uses the galaxy sequence and assumes that here the clus- 
ter members dominate, because we use the sequence galaxies 
within the central 400". In other words, we do not subtract the 
background and have enough statistics for a reasonable analy- 
sis. 

5. Summary and conclusions 

For the sample of clusters we discussed in this paper, the opti- 
cal luminosity we estimated, Tab. [5] correlates reasonably well 
with the luminosity observed in the X-ray band (data from 
Strubel & Rood, 1999). This is a well known result and we 
also observe that there is a tendency for the cluster luminosity 
to correlate with the redshift. Indeed by going at large distances 
any sample tend to pick up brighter objects. The cluster X-ray 
luminosity of the sample correlates as well with the redshift 
and that is also expected for the reason we just said. 



The important result we find, albeit with rather small sta- 
tistical significance due to the small span in luminosityes, is 
that the most luminous clusters in the X-ray band are compara- 
tively optically fainter than those having smaller luminosity in 
the X-ray band. The effect is rather robust if we consider the 
fact that it is manifest in all three fillers. Indeed in Fig. 1121 we 
plot the ratio L opt /Lx versus Lx and the effect is evident: the 
ratio L pt/Lx decreases as a function of Lx- The slope of the 
correlation, furthermore, decerases as a function of the wave- 
length, that is bright X-ray clusters are redder than fainter X-ray 
clusters. 



Fig. 12. The ratio of the optical Luminosity (in the g, r and i 
bands) of the clusters to their X-ray luminosity versus the X- 
ray luminosity. Brighter X-ray clusters are, percent wise, less 
luminous than fainter X-ray clusters. Furthermore the slope of 
the correlations seems to indicate that the cluster colors are a 
function of the X-ray luminosity in the sense that bright X-ray 
clusters are redder. 



The question arises whether the effect we observe is due to 
a transfer of matter between the galaxies and the intracluster 
medium (ICM) where the bright X-ray clusters had an higher 
stellar activity in the past. Some of these points will be dis- 
cussed in a subsequent paper. 

In summay, and over the years, we measured a set of clus- 
ters in various colors with the purpose of better understanding 
the shape of the luminosity function in clusters and to pin point 
the slope of the faint end. The observations were carried out 
in different colors and the analysis, and the interpretation of 
the data, made use of different algorithms to make sure that we 
were not biased by the method of analysis. All this work, while 
respectable, lead to a rather robust bi-modal cluster luminosity 
function, a result that is not new. Given the sample we used we 
were unable to estimate any variation of the LF as a function of 
the distance of the cluster center. While we find indication in 
the literature that there are variations (here we refer to the slope 
of the faint end), we feel that such results, since it is very impor- 
tant for the understanding of the cluster dynamics and cluster 
evolution, should be checked on large spectroscopic samples. 
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The important result of this work, however, consists in the 
estimate of the total luminosity in different optical color and its 
correlation with the X-ray luminosity. Indeed we believe that it 
is the estimate of the luminosity and colors of clusters that will 
give the opportunity, when compared to the X-ray luminosity, 
to better pin point some of the characteristics of the cluster and 
ICM evolution. 
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19.69!°* 


-1 34 +017 

-0 17 




i 


o.34:« : 5« 


16.69!°;^ 


0-88^:1? 


19.28!°| 


1 cn+0.07 
A -^^-0.07 




g 




i 7 99+O.24 
i / • ^-^-_q 24 


1.04^ 


in tq+0.30 

— 30 


_i 4 9+o.08 
-0 08 


A3667 


r 


16 +004 

-0.04 


16.66!™ 




18.67!°; 4 jj 


1 QQ+0.09 

—0.09 








16 32 +048 


1 03+ - 59 

-0.18 


18.22!™ 


1 ^9+0.07 
1 ■ J -0.07 




g 


0.8!°f 


18.52^| 


1 03 +023 

u -0.23 


19.69!°$ 


-1.63!°;! 4 


A3695 


r 


0-95^ 


17.90!°;! 




19.16!°;| 


-1-68!°;J° 




i 


0.38!°,$ 


17.88!™ 


0.90!°;| 


1 Q OQ+0.33 


_l C4+0.06 

w MJ.06 



Table 3. Parameters obtained with Maximum Likelihood method and the standard analysis 
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Cluster 


Filter 


G/S 


m G 


0~G 


m* 


a 




g 


14 +0 - 05 

-0 05 


17 14 +0 ' 40 

' -0 40 


i 79+0.31 

-0.31 


21 12 +a20 

-0.20 


_1 QQ+0.07 

— 07 


A4038 


r 


14+ 005 

-0.05 


16 54 +0 ' 58 


2.02!»;« 


20 - 88 -a23 


1 QQ+0.07 

-1.38_ 007 




i 


14 +0 - 04 

-0.04 


16.47!°;" 




20.86^!» 


_i 41+O.O6 

-0.06 




g 


1 16 +0 - 70 


18.20!^ 


1 -14+0.80 
-0 40 


19.06!|J0 


_j 7 4+o.i2 

' -0.12 


A4059 


r 


76 +4 2 ° 


17.22^ 


i oc + 1.40 

-0 60 


18.07+' °° 


-1 65 +a7 ° 




i 


i 7Q+4.20 
-0.30 


16 54+ 3 00 

-0.40 


i-os-oi!! 


17.87+J 2 ° 


-1 67 +on 

-0.11 




g 


55 +0 - 27 


17 - 82 -a? 6 


1.64! 2 6 f 


19.26^ 


-1 66 +006 


EXO 0422-086 


r 


0.43!»; 2 » 


17 - 46 ^24 


0.86!» 2 » 


18 - 98 -0 36 


_i 47+0.06 

-0 06 




i 


0.38^11 


17 26 +a28 

-0.28 


0.89^ 


18.39^ 


-1 40 +007 

-0.07 




g 


21 +0 - 85 


17 21 +0 ' 28 

-0.28 




18.97!^ 


-1 34 +017 

-0.17 


A0496 


r 


r\ 99+O.O6 


i6.so;:;: 


0.98*$ 


18.31^ 


1 AO+0.08 




i 


0-78^ 


16.46!» 2 | 


0.86!» : 2 » 


18.38!^ 


-1 4^+0.16 

-0.93 



Table 4. Parameters obtained with Maximum Likelihood method and the standard analysis 
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Cluster 


z 


Lx 


Filter 


Lo/L 


Lo/L 


Lo/L 


Lo/L 


L/Lx 


L/Lx 






lO^h^erg/s 




10 u rr 2 erg/s 

Flobs 


10 u rT 2 erg/s 

F lcom 


10 u rr 2 erg/s 
0.3Mpc 


10 u rT 2 erg/s 
A.R. 


Flobs 


A.R. 








g 


2.68 


2.81 


2.94 


3.68 


0.4954 


0.6802 


A0085 


0.0556 


5.41 


r 


2.11 


2.59 


2.69 


3.22 


0.3900 


0.5952 








i 


1.02 


1.69 


1.94 


1.94 


0.1885 


0.3586 








g 


2.49 


2.84 


2.86 


3.01 


1.1528 


1.3935 


A0133 


0.0566 


2.16 


r 


1.79 


1.92 


1.93 


2.03 


0.8287 


0.9398 








i 


0.74 


1.13 


1.14 


1.18 


0.326 


0.5463 








g 


3.46 


3.67 


3.72 


3.89 


0.7954 


0.8943 


A3667 


0.0556 


4.35 


r 


3.17 


3.25 


3.32 


3.58 


0.7287 


0.8230 








i 


1.96 


2.26 


2.27 


2.36 


0.4506 


0.5425 








g 


3.01 


5.63 


5.54 


5.97 


0.7099 


1.4080 


A3695 


0.0893 


4.24 


r 


2.44 


6.06 


5.97 


6.35 


0.5755 


1.4976 








i 


1.30 


3.86 


3.8 


4.04 


0.3066 


0.9528 








g 


2.08 


2.83 


2.84 


2.86 


2.0594 


2.8317 


A4038 


0.0292 


1.01 


r 


1.58 


2.93 


2.94 


2.95 


1.5644 


2.9208 








i 


1.05 


1.09 


1.09 


1.09 


1.0396 


1.0792 








g 


1.61 


1.64 


1.66 


1.72 


0.8994 


0.9609 


A4059 


0.0460 


1.79 


r 


1.59 


1.75 


1.76 


1.80 


0.8883 


1.0056 








i 


0.90 


1.03 


1.03 


1.03 


0.5028 


0.5754 








g 


1.58 


1.83 


1.86 


1.91 


1.2540 


1.5159 


EXO 0422-086 


0.0397 


1.26 


r 


1.19 


0.91 


0.92 


0.96 


0.9444 


0.7619 








i 


0.63 


0.85 


0.87 


0.92 


0.5000 


0.7302 








g 


1.37 


1.28 


1.37 


1.51 


0.6372 


0.7023 


A0496 


0.0328 


2.15 


r 


0.80 


0.99 


1.1 


1.25 


0.3721 


0.5814 








i 


0.37 


0.60 


0.61 


0.63 


0.1721 


0.2930 



Table 5. Luminosities obtained with the parameters of the standard analysis. X-ray luminosity computed in the 0.5 - 2.0 keV 
band assuming a power-law spectrum with energy index y = 0.4 (data from Strubel & Rood 1999). NOTE: the / luminosity 
of A0133, A3695, A4059 and EXO 0422-086 after background correction resulted to be slightly smaller than the luminosity 
of the few brightest galaxies, due to a sky background overcorrection and uncertainties in the cluster counts. In these cases, 
the luminosity was corrected assuming that brightest galaxies are cluster members as indicated by the morphology and by the 
redshift. 
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Abell 85 
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-1.5 -1.0 
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Fig. 13. Hereafter we show, for each cluster and for each filter used in this work, the fit of our model obtained with the maximum 
likelihood method and two boxes with the lcr, 2cr and 3cr confidence contours for the parameters of the Schechter (a vs. m*) and 
the Gaussian (mo vs. cr). In this page: A0085. 
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Fig. 14. A0133 
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Fig. 15. A3667 
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Abell 3695 Abell 3695 

3. or — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3. or 1 1 1 1 1 1 1 1 1 — 





Fig. 16. A3695 
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Fig. 17. A4038 
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Abell 4059 

3.0 f 1 1 1 1 1 1 1 1 1 1 1 r 



a. 5 - 
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Fig. 18. A4059 
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Fig. 19. EXO 0422-086 



